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ABSTRACT 

We report on our serendipitous discovery with the XMM-Newton Observatory of a luminous X-ray 
emitting cluster of galaxies that is located behind the Andromeda galaxy (M31). X-ray emission 
from the cluster was detected previously by ROSAT, and cataloged as RX J0046. 4-1-4204, but it was 
not recognized as a galaxy cluster. The much greater sensitivity of our XMM-Newton observations 
revealed diffuse x-ray emission that extends at least 5' and has a surface brightness profile that is well 
fit by the a-(3 model with (3 = 0.70 ± 0.08, a core radius = 56" ± 16", and a = 1.54 ± 0.25. A 
joint global spectral fit of the EPIC/MOSl, M0S2, and PN observations with Mewe-Kaastra-Liedahl 
plasma emission model gives a cluster temperature of 5.5 ± 0.5 keV. The observed spectra also show 
high significance iron emission lines that yield a measured cluster redshift of z = 0.290 with a 2% 
accuracy. For a cosmological model with Hq — 71 km s~^ Mpc~^, VLm = 0.3 and JIa = 0.7 we derive 
a bolometric luminosity of = (8.4 ± 0.5) x 10'''* erg/s. This discovery of a cluster behind M31 
demonstrates the utility of x-ray surveys for finding rich clusters of galaxies, even in directions of 
heavy optical extinction. 

Subject headings: galaxies: clusters: Intergalactic medium - X-rays: observation - Cosmology 



1. INTRODUCTION 

Galaxy clusters are the largest gravitationally bound 
structures in the universe. The evolution of cluster num- 
ber density of a given mass is sensitive to specific cos- 
mological scenarios (e.g. Press & Schechter (1974)). So 
observations of galaxy clusters are an important tool for 
constraining fundamental cosmological parameters. 

Due to the fact that 15% of the total cluster mass (e.g. 
Evrard (1997)) is in the form of hot diffuse plasma emit- 
ting at X-ray band via thermal bremsstrahlung (Sarazin 
1988), galaxy clusters are among the most luminous ob- 
jects in X-ray band. It makes X-ray selection an efficient 
means for constructing samples of galaxy clusters (see 
review by Rosati et al. (2002)) . X-ray selection has 
the advantage that the measurable X-ray luminosity and 
temperature are correlated with the cluster mass. 

Further, X-ray selection is useful for studying regions 
where optical searches are complicated because of dust 
extinction and heavy stellar confusion. X-rays are much 
less affected by extinction than optical photons and X- 
ray selection is almost free from source confusion prob- 
lems (Ebeling et al. 2002). Conducting X-ray selection 
based on i?OS'^T data at low Galactic latitude, Ebeling 
et al. (2002) were able to detect 137 galaxy clusters , 70% 
of which were new discoveries. 

With a new era of XMM-Newton and Chandra ob- 
servatories with their large effective areas and wide en- 
ergy ranges 0.3 — lOkeV, the capability of X-ray selec- 
tion increased. During an XMM-Newton observation of 
the galactic supernova remnant G21.5-09 located close 
to the Galactic Plane, Nevalainen et al. (2001) detected 
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a new galaxy cluster. Using only XMM-Newton data, 
they measured cluster redshift z — 0.1 to 1% precision 
that is especially important in regions with such strong 
optical source confusion, where the optical redshift mea- 
surements of galaxies are difficult . 

Here we present new XMM-Newton observations of the 
source RX J0046.4-|~4202 that indicate it is a high red- 
shift cluster located behind M31. RX J0046.4-I-4204 was 
detected during the first and the second deep ROSAT 
RSPC surveys of M31 performed in June 1991 and 
July/August 1992 respectively (Supper et al. 2001). 
Based on a comparison of the first and the second sur- 
veys, RX J0046. 4-1-4204 was classified as a potentially 
long term variable source. Our analysis of the data ob- 
tained with XMM-Newton revealed spatially extended 
emission, up to at least 5', from RX J0046. 4-1-4204. The 
observed spectra show iron emission lines that yield a 
measured redshift of z = 0.290 with a 2% accuracy. All 
these facts combined with optical image from Digitized 
Sky Survey allow us to conclude that RX J0046.4+4204 
is actually a distant galaxy cluster. 

In this paper, we assume the f2j\/ = 0.3 and = 
0.7 cosmology with the Hubble constant of Hq — 71 
km/s/Mpc. For the defined above cosmology and the 
measured redshift of z = 0.290, the angular size of 1' 
corresponds the physical size of 257 kpc. Statistical un- 
certainties are quoted at the 90% confidence level unless 
there is a statement saying otherwise. 

2. OBSERVATIONS AND DATA REDUCTION 

In the following analysis, we use the data from XMM- 
Newton observation of the XMM North 3 Field of M31 
centered at RA = 00''46™38^00; Decl = -h42°16'20.0". 
Fig. 1 (Left) shows optical image of M31 from Digitized 
Sky Survey with XMM-Newton FOV shown as a circle 
for M31 North 3 field. The XMM-Newton observation 
was performed on 2002 June 29 as a part of the Guar- 
anteed Time Program (PI: K.O. Mason). An analysis of 
the complete North 3 field is presented in Trudolyubov 
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Fig. 1.— Le/t.-Optical image of M31 from Digitized Sky Survey with JCMM-Newton FOV show as a circle for M31 North 3 field. Right: 
The combined MOS1-MOS2-PN vignetting-corrected image of the XMM North 3 Field of M31 in 0.8 - 2.5 koV energy band, square root 
intensity scale. 
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Fig. 2. — The Palomar Digitized Sky Survey optical image 
with overlaid brightness contours made from 0.8-2.5keV band 
MOS2/PN image. The X-ray image was smoothed with a 4" Gaus- 
sian kernel. The contour levels are based on the background noise 
going up in logarithmic steps. MOSl image was not used for con- 
struction of the contours because the center of the cluster falls at 
the edge of a CCD in MOSl camera. The angular size of the largest 
contour corresponds to ~ 0.7Mpc at 2 = 0.29. 

et al. (2005). 

In current analysis, we used the data from two EPIC- 
MOS detectors (Turner et al. 2001) and the EPIC-PN 
detector (Striider et al. 2001). The EPIC data was re- 
duced using the standard XMM-Newton Science Analysis 
System (SAS v6.1.0)^. We used the calibration database 
with all updates available prior to January, 2005. Only 
X-ray events corresponding to patterns 0-12 for MOS de- 
tectors and patterns 0-4 for PN detector were selected. 
All known bad pixels were excluded. 

^ See http://xmm.vilspa.esa.es/external/xmm_sw_cal/sas.shtml 



The EPIC background is highly variable and only its 
quiescent component can be accurately modeled. To de- 
tect and exclude the periods of high flaring background, 
we produced the light curves for each EPIC detector 
showing count rate in the 2.0-15.0keV energy band from 
the whole field of view but with detected sources masked. 
The light curves were binned to 200 s time resolution. We 
screened the EPIC data to recursively exclude time in- 
tervals with the deviation of the count rate exceeding a 
2(7 threshold from the average. Experiments with differ- 
ent choices of energy bands and flare detection thresh- 
olds have shown that our choice was close to optimal. 
The remaining good exposure time was ~ 24 ks for the 
EPIC-PN, - 41 ks for the EPIC-MOSl, and - 43 ks for 
the EPIC-M0S2. 

To account for strong XMM mirror vignetting, we used 
an approach proposed by Arnaud et al. (2001). Each 
photon was assigned a weight proportional to inverse vi- 
gnetting and these weights were then used in computing 
images and spectra. This was done using the SAS tool 
evigweight. 

Background modeling in our analysis was implemented 
following the double-subtraction method of Arnaud et al. 
(2002). The first step of this method is to subtract the 
particle-induced background component. This compo- 
nent can be estimated from a set of XMM observations 
with the filter wheel closed (so called "closed data"). 
We compiled the closed dataset from public observations 
available in the XMM data archive; these data were re- 
duced following identical steps as the science observa- 
tions. The closed background was adjusted to the cluster 
observations using the observed flux in the 10-15 keV 
band outside the field of view. The scaling factors are 
1.04, 1.05, and 1.01 for MOSl, M0S2, and PN respec- 
tively. The second step is to determine the cosmic X-ray 
background (CXB) component. Its spatial distribution 
should be fiat because vignetting correction is already ap- 
plied. Therefore, the CXB component can be measured 
in the source-free regions of the field of view at large radii 
from the cluster center (see § 3.1 and 3.2 below). 
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TABLE 1 

Results of Spatial Analysis. Combined 
EPIC-PN, MOSl AND MOS2 data, 0.8 - 2.5 
keV energy range. Parameter errors 
quoted are 90% confidence limits . 



Parameters 


P fit 


a -13 at 


a 




1.54 ±0.25 




0.60 ± 0.03 


0.70 ± 0.08 


rd") 


20.2 ±2.8 


56 ± 16 


rc(kpc for 2=0.290) 


87± 12 


240 ± 69 


xVdof 


156.7(119) 


139.4(118) 



Finally, wc applied a correction for photons registered 
during the CCD readouts, so called out-of-time events, 
to the EPIC-PN data. 

3. RESULTS 

The combined MOSl, M0S2, and PN image of the 
XMM North 3 Field of M31 in the 0.8 - 2.5 keV en- 
ergy band, corrected for the effects of instrumental vi- 
gnetting, is shown in Figure 1 (Right). The raw image 
was convolved with a Gaussian function with spatial scale 
fj = 4". We define the cluster center to be at the location 
of the X-ray surface brightness peak of the cluster emis- 
sion in the 0.8 - 2.5 keV energy band, a = 00''46'"24^8 
S = +42°04'26" (J 2000), with an estimated uncertainty 
radius of 6" (90% CL), determined by the wavelet de- 
composition algorithm of Vikhlinin et al. (1998). The 
Palomar Digitized Sky Survey optical image shows no 
extended optical counterpart for RX J0046. 4+4204 (sec 
Fig. 2). 

During the following spectral and spatial analysis, we 
excluded all detectable point sources from the data. The 

sources were detected separately in the "optimal" 0.3- 
3 keV, "soft" 0.3-0.8 keV, and "hard" 2.0-6.0 kcV energy 
bands using the task wvdecomp of the ZHTOOLS package 
^. Detected point sources were masked with circles of 
80% PSF power radii. 

3.1. Spatial Analysis 

Spatial analysis of the cluster emission was performed 
in the 0.8 — 2.5 keV energy band using all detectors. Ex- 
periments with different energy bands showed that the 
signal to noise ratio was close to optimal for the chosen 
band. We used a pixel size of 4" in our spatial analysis. 
The image for each camera was corrected for vignetting. 
The PN image was corrected for out-of-time events. We 
subtracted the particle background component from the 
images as described in § 2. 

We extracted the azimuthally averaged surface bright- 
ness profiles centered on the X-ray surface brightness 
peak, excluding the CCD gaps and circles around the 
point sources. The profiles were logarithmically binned 
with a step of Ar = O.lr. A logarithmic radial bin- 
ning approximately preserves the signal to noise ratio in 
annuli until the background becomes comparable to the 
signal. For the extracted profiles, the chosen step keeps 
the signal to noise ratio above 3 in annuli until ~ 5'. The 
obtained profiles were used to derive the parameters of 

® Seehttp: / /hea-www.harvard.edu/ ~alexey/zhtools/ 



the spatial distribution of the ICM, cluster fluxes, and 
the level of the CXB component. 

The cluster surface brightness profiles are often mod- 
eled with the so called /3-model, oc (1 -|- r'^/r1)~^^ 
or Sx oc (1 -)- 7.2/y.2-)-3/3+o.5 (CavaHere & Fusco-Femiano 
1976). However, this model poorly describes clusters 
with sharply peaked surface bigness profiles related to 
the radiative cooling of the ICM in the cluster centers. 

We used the simple modification of the /3-model 



" (1 + rVr2)3/3-"/2 ' ^ ^ 

which allows to model a power law-type emission excess 
in the cluster centers. The similar models were also used 
by Pratt & Arnaud (2002) and VikhUnin et al. (2005). 

The model for the observed surface brightness profiles 
can be obtained by numerical integration of eq.(l) along 
the line of sight, and convolution of the result with the 
XMM PSF^. To represent the uniform sky X-ray back- 
ground, we added a constant component to the model 
and treated it as a free parameter. The values of a, (3, 
and Tc were derived from the joint fit to the observed 
profiles in the MOSl, 2 and PN cameras, with the overall 
normalizations and background levels fit independently 
for each camera. The obtained parameters for the a-/3 
model are summarized in Table 1. 

We obtained a (3 value of 0.70 ± 0.08. For a sample of 
local clusters, Vikhlinin et al. (1999) showed that if cool- 
ing flow regions were excluded from /3-model flt, values 
of the [3 parameter were distributed over a narrow range 
0.7 ± 0.1. So the value of the f3 parameter we derived 
describes the typical distribution of the ICM. 

The obtained value of the core radius , rc = (56 ± 16)" 
((240 ± 69) kpc for z==0.290 ), is similar to the values of 
the Tc parameter Vikhlinin et al. (1999) measured for a 
sample of local clusters. 

For the value of a parameter we obtained a = 1.54 ± 
0.25. This value is consistent with values of the a param- 
eter Vikhlinin et al. (2005) derived for typical cooling- 
flow clusters using a sample of local clusters. Vikhlinin et 
al. (2005) in their analysis used Chandra data where the 
PSF is not an issue, so their measurements of a param- 
eter are more direct and do not depend on the quality of 
a PSF model as in our case. 

While the EPIC-pn PSF is azimuthally symmetric and 
is well calibrated up to large offset angles in the 0.8-2.5 
energy band (see Fig 1-6 in Kirsch (2005)), the EPIC- 
MOS PSF is less reliable and seems not azimuthally sym- 
metric. The asymmetry is not currently modeled in the 
SAS calibration. To check how imperfect calibration of 
the EPIC-MOS PSF can influence the obtained result, 
we fitted the surface brightness profile of the cluster ex- 
tracted from the EPIC-PN detector alone. The obtained 
best-fit values were within our statistical uncertainties. 

We can quantify the strength of cooling flow calcu- 
lating the fraction , /yokpc of the count rate coming 
from a central region of fOkpc to the total count rate 
(Peres et al. 1998). This quantity exhibits large disper- 
sion, but in general, clusters with massive cooling flows 
show /rokpc > 15% (Peres et al. 1998). For a typi- 
cal non-cooling-flows cluster /yokpc ^ 6% (Markevitch 

^ Wc used the latest available values of the King function 
parametrization of the MOSl, MOS2, PN PSF, see 
http://xmm.vilspa.esa.es/docs/documents/CAL-TN-0018-2-4.pdf 
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1998). Wc estimated the count rate coming from the 
70kpc central region using the best fit a-f3 modcL The 
total count rate was calculated integrating the observed 
surface brightness profile within 1.4Mpc radius. We ob- 
tained /rokpc = 18%. 

The combined M0S1-M0S2-PN surface brightness 
profile along with the best fit a-f3 model are shown in 
Fig. 3. From Fig. 3 one can see how the XMM PSF can 
flatten a peaked profile. For comparison, we also fit the 
surface brightness profiles by the standard /3-model set- 
ting a = 0. The obtained parameters are summarized in 
Table 1. 

3.2. Global Spectrum 

For our spectral analysis, we used the screened 
vignetting-corrected data in the 0.5 — 10.0 keV energy 
band from all cameras. The spectra of the cluster were 
extracted from a circular region with angular radius of 
3.1' (O.SMpc for 2=0.290) for aU EPIC data. AU de- 
tected contaminating point-like sources were excluded 
from the source and background extraction regions. We 
subtracted the particle background component from the 
images as described in § 2. To estimate the remaining 
CXB component, we extracted the spectra from a cir- 
cular region with angular radius of 9.0' centered at the 
on-axis position , but a part of the region falling into the 
9.65' (2.5Mpc for z=0.290) circle centered at the cluster 
center was excluded. The response matrices and efli'ective 
area flies were generated by the standard SAS tasks. Be- 
cause the data were previously vignetting corrected, the 
effective area files were created for the on-axis position 
using the routine arfgen. The response matrices were 
generated in the spectrum extraction region via rmfgen. 

The source spectra were binned to have at least 30 
counts in each spectral bin and fit in XSPEC 11.3.0 (Ar- 
naud 1996) by the Mewe-Kaastra-Licdahl plasma emis- 
sion model (Mewe et al. 1985). We used abundances 
from Anders & Grevesse (1989). Galactic photoelectric 
absorption was accounted for using the WABS model 
(Morrison & McCammon 1983). The spectra from the 
EPIC-PN (3677 counts), MOSl (1666 counts) and M0S2 
(2179 counts) detectors were fitted both jointly and sep- 
arately. For the joint fits, only spectral model normaliza- 
tions were allowed to vary independently. The results of 
both joint and separate spectral fitting of the EPIC-PN, 
MOSl and M0S2 data are summarized in Table 2. 

We obtained a redshift value of 0.290 with 2% accuracy. 
The redshift values estimated independently from EPIC- 
PN, MOSl and M0S2 data are in good agreement within 
measurement errors. For the value of hydrogen column 
density we obtained Nh = (2.2±g;^) x 10^^ cm that is 
significantly above the Galactic hydrogen column in the 
direction of M31, ~ 7 x 10^°cm~^ (Dickey & Lockman 
1990). We checked that the values of absorption obtained 
independently from EPIC-PN, MOSl and M0S2 data 
are consistent within the measurement errors. 

The EPIC-PN, MOSl and M0S2 spectra, along with 
the best-fit spectral models, are shown in Figure 4. Note 
that separate spectral fitting of the EPIC-PN, MOSl and 
M0S2 data gives consistent values of the model param- 
eters. 

3.3. Temperature profile 



To construct temperature profile, we extracted indi- 
vidual spectra in 5 annuli: — 50kpc , 50 — lOOkpc , 
100 - 200kpc, 200 - 400kpc, and 800 - 1600kpc, centered 
on the position of the X-ray surface brightness peak us- 
ing data from all cameras. After background subtraction, 
the number of counts (M0Sl-hM0S2-hPN) in 0.5-lOkev 
energy band in 1, 2, 3, 4, and 5 annuli was 632, 1275, 
2071, 2069, 1063, and 459 accordingly. 

The correction for the XMM PSF effect was done us- 
ing an approach of Pointecouteau et al. (2004). Using 
the best-fit a-/3 models of the cluster brightness and the 
XMM PSF calibration, we calculated the redistribution 
matrix, Rij ,of each temperature to each annulus which 
represents relative contribution of emission from annulus 
i to the observed flux in annulus j. The model spectrum, 
Sj, is then given by 

5, =^i?,,^(T,), (2) 

where Ti is the temperature in annulus i and S{Ti) is 
the mekal spectrum for this temperature. Fitting this 
model to the observed spectra in all annuli simultane- 
ously and treating all Ti as free parameters gives the de- 
convolved temperature profile. Unfortunately, the statis- 
tic was poor to measure metallicity profile. So for all 
annuli we fixed the metallicity values at the best-fit val- 
ues obtained from the global spectrum. We checked that 
allowing the absorption and the metallicity to be freely 
fitted does not significantly change the result. The val- 
ues of redshift and absorption values were also fixed at 
the best-fit values. 

We fitted the observed temperature profiles by a 3-D 
temperature model : 

T(r) = To/(l. + {x/0.6fr, (3) 

where x = r/rg. A similar model describes the tem- 
perature profiles for low redshift clusters (Vikhlinin et 
al. 2005). For local clusters, ro scales with the clus- 
ter temperature as tq ^ 0.50 T^/^ Mpc where T is in 
keV (see Fig. 16 in Vikhlinin et al. (2005)). For 
the best-fit temperature of T = 5.5keV derived from 
the global spectrum, ro = 1.17keV. We fixed ro at 
the value suggested by the low redshift clusters only 
with an additional redshift scaling ro oc l/E(z), where 

E(z) = y^^Aiiz + l)^ + ^A, to account for the redshift 
dependence of the virial radius for clusters with a flxed 
temperature (Bryan & Norman 1998). To fit the ob- 
served profiles, we projected the 3-D model along the 
line; of sight using the emission measure profile from the 
best flt a-P model. Projection was based on a weighting 
method which correctly predicts the best-fit spectral T 
for a mixture of different temperature components (Maz- 
zotta et al. 2004; Vikhhnin 2005). 

Our 3-D model has a flat shape in the central region 
(r < 0.6ro/7°'^) and is not intended for modeling central 
temperature drops associated with a cooling region. The 
innermost bin shows an apparent temperature drop and 
can indicate the presence of a cooling flow. Therefore we 
excluded the central bin from the fit. Taking into account 
the fact that the temperature profile was corrected for the 
XMM PSF this procedure is correct. The obtained best- 
fit model, the deconvolved and original temperature pro- 
files are shown in Fig. 5. The uncertainties on the best- 
fit model were calculated from Monte-Carlo simulations. 
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Fig. 3.— Left: Upper panel: The combined MOS1-MOS2-PN surface brightness profile of RX J0046.4+4204 in the 0.8 - 2.5 keV energy 
band. Black solid line shows ct-fi model best fit convolved with the XMM PSF. Black dashed line shows o-/3 model best fit. Black dotted 
line shows the XMM PSF. Black dotted-dashed line shows the level of the cosmic X-ray background component. Lower panel: The residual 
between the data and the best-fit model in terms of sigmas. Right: The same with the standard /3-modcl best fit. 
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Fig. 4.— Left: Count spectra of RX J0046.4+4204 taken with the XMM-Newton/EPIC-PN (black circles), MOSl (blue triangles) 
and MOS2 detectors (magenta squares). The corresponding best-fit spectral models (absorbed red-shifted Mekal model) convolved with 
instrumental responses are shown as black, red andblue solid lines. Right: Expanded view of the 4.0-7.0keV energy band. Upper, Middle 
and Lower panels show PN, MOSl and MOS2 data. A red-shifted iron emission line feature is clearly evident. 



We applied Gaussian scatter to the observed tempera- 
ture profile within its uncertainties and fitted simulated 
profiles. The best-fit model uncertainties estimated as 
rms scatters in narrow radial bins from 1000 simulations 
are shown in Fig. 5.. 

3.4. Mass and Luminosity measurements 



Assuming hydrostatic equilibrium for the ICM, we can 
use the best fit temperature and density profiles to derive 
the total cluster masses: 

Mf ^ rr(r) f dlogpjr) ^ d\ogT{r) \ 



GfiTTip \ dlogr 
We calculated the total mass of M^qq 



d log r ^ 

(5.1 ± 1.1) X 
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TABLE 2 

Results of Spectral Analysis. EPIC-PN, MOSl and MOS2 data, 
0.5 — 10.0 keV energy range. Spectral extraction radius is 3.1'. 
Parameter errors quoted are 90% confidence limits . 



Parameters 


MOSl 


MOS2 


PN 




Combined 


kT(keV) 

Zq 
z 


0.21±^| 
6.4±p 
0.44±i:i 
0.296±«:W^ 


72-1-8; 5o 


0.21±0 
5.1±» 
0.45±S 
0.287±[] 


02 
92 

4 

008 


n 22+''.02 
B.02 

e K-l-0.5 

Q-5 

0.57it 
290+^'-^^5 




0.01 0.1 1 

R,Mpc 



Fig. 5. — Upper panel: Same as Fig. 3 (Left upper panel). Lower 
panel: The temperature profile of the cluster as a function of an- 
gular radius obtained from EPIC-PN/MOS1/MOS2 data. SoUd 
black circles show the deconvolved projected temperature profile. 
For comparison, open circles shows the raw measurement from the 
X-ray fit in the same annuli.The error bars correspond to 68% (l-c) 
confidence limits. Solid line shows the best-fit projected tempera- 
ture profile and dashed lines correspond to its 68% CL uncertain- 
ties. 

1O"M0 and the radius of rgoo = 1.09 ± O.OSMpc corre- 
sponding to the mean overdensity A — 500 relative to 
the critical density at the cluster redshift. The uncer- 
tainties on the total mass due to the temperature profile 
were calculated analogically to the uncertainties on the 
best-fit temperature profile model. We calculated the 
total mass for each simulated profile. Then the uncer- 
tainties were calculated as the boundaries of the region 
containing 90% of all realizations. The total mass un- 
certainties due to the error on the density gradient were 
calculated following Pratt & Arnaud (2002). The value 
of c? log p(r)/(i log r at rsoo was considered as an indepen- 
dent parameter of the a-f3 model instead of /3. We refit 
the surface brightness profile with the new parameter 
set and measured the uncertainties on d log p(r)/(i log r 
at ^500. The final total mass uncertainties were calcu- 
lated by adding quadratically the total mass error due 
to the temperature profile and the density profile. The 
uncertainties on i?5oo are related to the total mass un- 
certainties as CTMboo 7*^500 = ScTrsoo 7^500 ■ 

We calculated the emission measure- weighted temper- 
atures (volume-averaged with weight w — p^) of Temw = 
6.0 ± 0.7 within 7Gkpc < r < rsoo- It is not a surprise 



that obtained Tgmw is higher than the best-fit spectro- 
scopic temperature because we excluded the innermost 
bin from the temperature profile fit. The uncertainties 
for Temw were calculated from Monte-Carlo simulations 
in the similar way as for the temperature profile best-fit. 

We can compare the obtained M500 with predicted 
value from M — T relations. We will use as a low- 
redshift reference the M — T relation from Vikhlinin et al. 
(2005), which is similar to the M — T relation measured 
by Arnaud et al. (2005). The M -T relation derived by 
Vikhlinin et al. (2005) predicts M500 = 6.1 x lO^-^Mo for 
a 6 keV cluster. The self-similar theory predicts that for 
the same temperature the mass evolves as A^a oc E{z)^^ 
with z, where E{z) = H{z)/Hq = (0.3(1 + z)^ + 0.7)^/^ 
(e.g. Bryan & Norman (1998)). Dividing the predicted 
mass by E{z = 0.29) to place it at z = 0.29, we obtain 
-^500 = 5.3 X This value is in agreement with 

the derived M500 = (5.1 ± 1.0) x W^^Mq. 

We derived the unabsorbed bolometric luminosity of 
= (8.4 ± 0.5) X 10'^'^ h^^ erg s-^ using the Mewe- 
Kaastra-Liedahl plasma emission model. We used Temw 
as the temperature parameter for the model. The model 
was normalized by the following count rate: we sub- 
tracted from the observed 0.8 — 2.5keV count rate cal- 
culated within < r < 1.4Mpc the count rate calculated 
from the best-fit a-/3 model within r < TOkpc, and mul- 
tiplied the result by 1.06 (Markevitch 1998). 

The fact that our analysis is similar to the one done 
by Markevitch (1998) for low redshift clusters allows us 
to compare the obtained luminosity with the prediction 
from their L — T relation. The L ~ T relation derived 
by Markevitch (1998) predicts the luminosity of 6.2 x 
10'^'' h^^ erg for a 6 keV cluster. It is well below the 
derived luminosity. However, correcting the predicted 
luminosity for the evolution in L — T relation, = {1 + 
z^-^Lq, reported by Vikhlinin et al. (2002), we obtain 
9.2 x lO^** h^i erg s~^. This value is close to the observed 
one. 

4. DISCUSSION AND CONCLUSIONS 

Our deep XMM-Newton observations of M31 have 
shown that RX J0046. 4-1-4204 is not located in that 
galaxy, but rather is actually a distant cluster of galax- 
ies. We found RX J0046. 4-1-4204 has spatially extended 
X-ray emission and that the spectrum clearly shows a 
red-shifted iron emission line. Straightforward fitting of 
the iron line yields a cluster redshift of Z=0.290 with 
2% accuracy and that the redshift values estimated in- 
dependently from EPIC-PN and MOS data are in good 
agreement within measurement errors. 

The large scale spatial distribution of RX 
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J0046.4+4204 is well fit by the model with 

P = 0.70 ± 0.8, a core radius = (56 ± 16)" or 
rc = (240 ± 69)kpc for z=0.290 , and a = 1.54 ± 0.25. 
The obtained values of f3 and rc are consistent with the 
parameters of typical clusters (Vikhlinin et al. 1999). 
The derived a agrees with values measured for cooling 
flow clusters (Vikhlinin ct al. 2005). 

The spatially integrated X-ray continuum is well fit 
by red-shifted (z=0.290) Mewe-Kaastra-Liedahl plasma 
emission model with low energy photo-electric absorp- 
tion. The best fitting global model for the joint PN, 
MOSl, and M0S2 measurements yields the parame- 
ters kT{keV) = 5.5 ± 0.5, fractional solar abundance 
Zq = 0.57±°a:\l, ,a redshift z = 0.290 ± 0.005, and a 
column depth = 2.2(±q;2) x lO^^cm"^. This derived 
column depth is significantly larger than the Galactic hy- 
drogen column in the direction of M31, ~ 7 x 10^*^ cm~^ 
(Dickey & Lockman 1990). 

It is interesting to compare the value of absorption 
obtained for the cluster with absorption for the nearby 
X-ray sources. Fig 1 (Right) shows five bright X-ray 
sources lying in the direct vicinity of RXJ0046. 4+4204 
for which Trudolyubov et al. (2005) were able to measure 
column depth based on their X-ray spectra and to iden- 
tify some of them. XMMU J004540.5-h420806(source 
#5 in Fig. 1) was identified as a foreground star 
with the column depth of 0.2 x lO^-'^cm"^ (two sigma 
upper hmit). XMMU J004648. 0+420851 (source #3) 
was identified as a background radio source with 
the column depth of (4.0±^g) x lO^icm-^. XMMU 
J004627.0+420151 (source #1) was identified as a glob- 
ular cluster source in M31. The derived column depth 
for this source was (1.3 ± 0.1) x lO^^cm^^. The nature 
of two last sources, XMMU J004611.5+420826(source 
#4) and XMMU J004703.6+420449 (source #2), was 
unclear. Trudolyubov et al. (2005) proposed that 
these sources could be two AGN located in the back- 
ground of M31. The derived column depths of XMMU 
J004611. 5+420826 and XMMU J004703.6+420449 were 
(2.5 ± 0.3) X 102icTO-2 and (2.3 ± 0.6) x lO^icm-^, re- 
spectively. 

The absorption value of the background radio source 
is 2 times higher than the cluster absorption and most 
likely to be intrinsic. The column depth of the glob- 
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ular cluster candidate is smaller than that of the clus- 
ter. To explain this fact, it might be proposed that the 
globular cluster candidate is located in front of the disk 
of M31, while the cluster is obscured by the disk. On 
the other hand, XMMU J004611.5+420826 and XMMU 
J004703. 6+420449 have column depths consistent with 
that of the cluster, suggesting that they could be also 
located in the background of M31, although it is unclear 
what fraction of the column depths is intrinsic. This in- 
terpretation is in general agreement with results of spec- 
tral fitting of a large sample of M31 globular cluster X- 
ray sources (Trudolyubov & Priedhorsky 2004). Tru- 
dolyubov & Priedhorsky (2004) found that globular clus- 
ter sources located in front of M31 disk have typical val- 
ues of absorbing column in the range of (0.5—1.5) x 10^^, 
while the sources located behind the disk or embedded 
into it show higher absorbing columns of (2 — 4) x 10^^ 
cm 

The extracted temperature profile corrected for the 
XMM PSF shows the central temperature decline that 
confirms the indication of cooling fiow presence from 
the spatial analysis. Using the spatially resolve tem- 
perature profile we derived -R500 = 1.09 ± 0.07Mpc, 
M500 = (5.3+1.0) x10"Mq, andTe^^ = 6.0+0.7 within 
i?50o- The values of M500 corrected for the evolution and 
Temw are in agreement with local M — T relations. 

The study we have presented here shows the utility of 
sensitive X-ray observations for identifying and studying 
clusters of galaxies in directions where foreground confu- 
sion or heavy optical extinction makes optical selection 
complicated. 
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